Stability of quantum breathers 
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Using two methods we show that a quantized discrete breather in a 1-D lattice is stable. One 
method uses path integrals and compares correlations for a (linear) local mode with those of the 
quantum breather. The other takes a local mode as the zeroth order system relative to which 
numerical, cutoff-insensitive diagonalization of the Hamiltonian is performed. 
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Discrete breathers pj are nondispersive classical exci- 
tations that are known to be significant in a variety of 
natural systems. Recent reviews 0,01 give examples and 
references. The quantum theory of (discrete) breathers 
does not yield easily to numerical simulation and is not 
as fully developed. 

Over the past decade there has developed a body of 
literature £| in which it is claimed that unlike classical 
breathers, the quantum version cannot in principle be 
stable. This conclusion is reached by considering the ef- 
fect on the phonon field of the rapidly oscillating particle 
at the core of the breather. Lifetimes have been calcu- 
lated, and for alkali halides are about 10 ns. Other work, 
also treating the breather as an external classical force, 
0, challenges this estimate and suggests stability in the 
limit of large systems. 

A fully quantum result 6] found eigenstates of a model 
Hamiltonian with a dropoff in correlation function that 
is strongly suggestive of breathers, although translational 
invariance prevents complete localization. As recognized 
in [f| (and below), such approaches require cutoffs which 
may affect the conclusion. 

In this article we make a two-pronged attack. We do 
a numerical diagonalization, similar to that of |6|, incor- 
porating methods to reduce cutoff effects and to provide 
direct evidence of localization. In addition we do a path 
integral calculation, treating the breather fully quantum 
mechanically, similar to Feynman's polaron The ap- 
proximations are different: in one case a cutoff, in the 
other semiclassical asymptotics. 

The Hamiltonian in is 



±LJ 2 T 2 
2 UJ s x k 



+ \^x)l + \\xt} , (1) 



with xq = xn+i, (Ax)fe = Xk — Xk+i, and Xk G K. Wc 
also consider interactions of the form ^2{&-x)\, which 



resembles our own model J8J. An important simpli- 
fication is to remove all nonlinear interactions except 
those of particle-0, the one with large amplitude motion. 
This strategy has been adopted by other authors and 
we checked whether this affects the general properties of 
the classical breather. It does not a consequence of 
the remarkable breather property, namely that a central 
atom vibrates strongly, while its neighbors hardly move 
(so that nonlinear forces are negligible for them) . Claims 
that the breather decays are based on this same approx- 
imation. 

We emphasize the reasoning leading from a transla- 
tionally invariant breather to what is effectively a non- 
linear local mode. Although quantum tunneling requires 
that the breathers form a band, information on the clas- 
sical breather corresponding to our excitation — including 
in the presence of moderate noise — indicates that tunnel- 
ing is small, allowing us to drop translational invariance. 
This perspective has been adopted by others as well 0,0- 
Having done this, we go beyond previous work by quan- 
tizing the breather atom itself. In other approaches one 
has de facto a local mode, interacting with the lattice. 
We retain the nonlinear interactions (implicitly), while 
also quantizing the principal breather atom. 

In both our approaches we make use of a linear local 
mode. If in (the modified) Eq. , one replaces Xxq /4 by 
lu 2 Xq/2 with large the resulting system is well under- 
stood both classically and quantum mechanically. You 
have both classical confinement and quantum stability. 
Our use of the local mode parallels others' treatment of 
the quantum breather as a phonon field in the presence 
of rapid classical oscillations. The impact of the breather 
cannot be considered small. They treat that impact by 
means of a classical oscillating field. We use the quantum 
local mode. For the path integral, that use is the com- 
parison with local mode correlation functions. For the 
diagonalization scheme we perturb around that mode. 
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Path integral approach. The path integral allows the 
elimination of quadratic degrees of freedom at the ex- 
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pense of introducing nonlocal self coupling. Dropping 
quartic terms for all but Xq, the Lagrangian is 



C = E,to{in-^o(Ax)„-^}/2-A4/4 
-fix (x 



(2) 



The fictitious coupling, nxo(x m + Xff+i— m ), allows study 
of localization for m far from 0. The derivative (at /i = 
0) of the propagator provides a (O-m)-correlation. By 
standard methods 0, El El we path integrate degrees 
of freedom 1—N, which form a chain with known normal 
modes. These modes see a forcing term from the so-far- 
unintegrated Xo(t). The propagator, G, is a function of 
the endpoints of all Xk- One takes the matrix element 
of G in the chain ground state (a la Q) and divides by 
corresponding free chain matrix elements. The result is 



G(q f ,T;q 



<,o) = / 



Vqe 



*(So+S cff ) 



(3) 



where q = xq, the action Sq arises from the original C 
sans chain terms, and S et{ comes from the integration 
over chain degrees of freedom. Specifically 



So = J dt {q 2 /2 - (cj 2 /2 + ul) q 2 - Xq 4 /l} 



(4) 



= / / dtdsK(\t-s\)q(t)q(s), (5) 
o Jo 



with 



T 2 = 



iV + 1 



oj n Sin ■ 
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/i sin 



AT + 1 



(6) 
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and f2 2 = a; 2 + 4w 2 sin 2 [mr/2 (N + 1)] is the spectrum of 
the chain. 

We study Q in the stationary phase approximation, 
that is we evaluate So + S Bfi along extremal "classi- 
cal paths." In principle one can get the ground state 
energy by going to large imaginary times, for which 
Q(q, —iT; q) ~ |<M<?)| 2 ex P( — TEq). In practice S — > const 
and d 2 S/dq l dqf gives this exponentially small quantity 
(as for the harmonic oscillator). Such precision was not 
numerically possible; in any case, knowing the energy 
does not establish localization. Note that <f>(q) is the 
overlap of the ground state of the (N + l)-atom ring with 
that of the A^-atom chain. 

To study localization we compare Q with the propaga- 
tor of a corresponding local mode, i.e., a system in which 
Xq A /A is replaced by uj\q 2 /2 for appropriate u\ (~ vA). 
To derive Qi oaal we redo the process described above. K is 
the same. For this system we know that as u>i increases, 
the vibrations are localized near 0. How can that be seen 



in Siocai or S7 The method is to go to imaginary time 
and vary /j, (the coupling to atom #to w AT/3). 

Before presenting results we comment on numerics. 
The extremal solves a two-time boundary value prob- 
lem with a nonlocal interaction. In the linear case we 
discretize time, so that d 2 /dt 2 becomes a difference ma- 
trix, nearly diagonal, but the AT-matrix is spread out. 
Call the resulting operator B = [d 2 /dt 2 ] + .... Take 
q to be a column vector, qj <-> q(tj). For q(0) = a, 
<7(T) = b, and e = tk+i — tk, the extremal is q = B~ 



where qf ] = -(a, 0, . . . , 0, b)/e 2 (cf. [13). Note that 
S = qq/2\J l . To deal with computer limits on e, we calcu- 
lated S for the smallest practical e and for larger values. 
S as a function of e was then extrapolated to e = 0. 
This technique was validated for the demanding energy 
calculation, involving the exponentially small d 2 S/dadb. 
Here the answer is known from the classical chain and 
ring frequencies. (Eq also requires a T — > oo extrapola- 
tion.) The function 0, a Gaussian whose spread can be 
calculated by standard methods, provided another check. 
Numerically that spread can be extracted by calculating 
S with varying endpoints and fitting the slope of \/S. 
This too confirmed the method. 

With q* the nonlinear, nonlocal two-time boundary 
value problem cannot be solved by matrix inversion. Our 
method was a variation of , defining a subsidiary vari- 
able z(t) = J Q Ki(\t — s\)q(s)ds, where the subscript 
on K\ indicates that K (= iK(-iu) with T — > — iT) 
is replaced by a single "cosh" with optimized parame- 
ters, in fact an excellent approximation. This allowed 
a local boundary-value numerical technique, with a self- 
consistency demand on z(0). Because of the approxima- 
tion we further attempted small variations to lower the 
true action, but they resulted in essentially no change. 

Fig. H shows the results of changing n for the local 
mode (a) and for the breather (b). Localization is de- 
duced from the variation of dS/dfj, as a function of both 
a {= q' = q 1 ) and the coupling, u>\ or A. All runs are for 
T = 2, which (given the energy scale) allows many states 
besides the ground state to survive in a spectral sum for 
Q. Therefore a large a boundary condition selects for the 
excited states, in particular for a localized excited state 
(breather or local mode) if there is one. When the cou- 
pling is small, there is no localized state; forcing q to be 
large forces all atoms on the ring to depart from their 
usual positions. Hence the large magnitude of dS/d^i on 
the left of both figures — for the largest a values. (Small 
a has almost no effect even for small coupling.) For large 
coupling however, there is a localized state: both for the 
linear (u>\) and nonlinear (X) systems, forcing q to depart 
from has almost no effect on a distant atom. As fur- 
ther checks, in we show that dS/dfj, — + for A — > oo, 
and that for both linear and nonlinear cases dS/d/i has 
essentially the same behavior as T grows El . 

Diagonalizing a truncated Hamiltonian. The Hamilto- 
nian Q is made finite-dimensional by using a phonon 
basis and imposing a cutoff on the level of phonon exci- 
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FIG. 1: dS/d/j, (essentially a correlation function) has 
markedly different behavior for large and small endpoint ("a") 
values for q. Both for a linear local mode (as a function of 
u>i) and for the breather (function of A), a large value of a de- 
mands relatively large correlations when neither local mode 
nor breather is present (uii ~ A ~ 0), but that correlation is 
wiped out for large parameter, for which the breather or local 
mode is effectively decoupled from the rest of the ring. 



tation. A low cutoff is needed because of the proliferation 
of dimensions in the (implicit) tensor product of phonon 
operators. The same problem was faced in [f| and we al- 
leviate it by the following strategy. Instead of perturbing 
around free modes, we use modes of a fictitious (linear) 
system with a local mode at site-0. This reduces the am- 
plitudes of other modes in the quantum breather. Eq. IjTjI. 
with nonlincarity confined to xq, is written 



H 



N 

£ 

k=Q 



{H + Ws*l + ^U^)l}+^i4+vi , (8) 



with Vi = Axq/4 — loIxq/2 the perturbation. A varia- 
tion on Eq. Q uses a nonlinear coupling, that is, terms 
\{Ax)\/A. Here too, we drop nonlinear terms except 
those (two) that include xq, and add and subtract corre- 
sponding quadratics. As for the classical dynamics 
is substantially the same. 

We outline the calculation: The classical local mode 



TABLE I: Principal components of the breather state. The 
first 4 columns refer to the 4 symmetric phonons in a 6-atom 
ring. (N atoms => [N/2] + 1 symmetric modes.) Row 1: fre- 
quencies (local mode is highest). Subsequent rows: number- 
operator values. Fifth column: norm squared of the mode 
(log 10 in parentheses). Cutoffs: local mode 13, others 6. For 
this state first order perturbation theory is good to 0.3%. The 
last column reports the same calculation with cutoff 8. 
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problem is solved, giving rise to (truncated) creation and 
annihilation operators in terms of which the perturbation 
is expressed. The operator xq — XX a ^ + a ^) u ^(0) / 'y / 2H,£, 
where ue is the ^-th mode of the system and I = is the 
local mode. The key to wi's effectiveness lies in {up(0)}. 
For io\ = these are all of about the same magnitude (for 
N = 20 they are ~ 0.3 for those not zero by symmetry). 
But with u>i — 3, u (0) ~ 0.988, while the (nonzero) oth- 
ers average ~ 0.04. This severely reduces the amplitudes 
of phonons other than the local mode in the eigenstates 
of the full Hamiltonian. 

With A = 8, Wo = 1) and uj s = 1 (used in @), and for 
lui = 2.5 the true eigenfunction having greatest overlap 
with the first excited state of the local mode is shown 
in Table IB Clearly the local mode dominates. The next 
largest component is the thrice excited local mode, which 
is merely a shape adjustment. Other modes barely make 
the 10 -3 level. Note that the highest excitation level for 
other phonons is 3, indicating that a cutoff of 6 is safe. 
In fact, even to probability 10 -8 there is no excitation 
higher than 3 except for the local mode. 

As a check of cutoff sensitivity we repeated this calcu- 
lation with a cutoff of 8 (but the local mode still at 13). 
The results are in the last column of Table HJ There is 
little sensitivity to the change — not only in the probabil- 
ities but in the composition of the state. Increasing ring 
size and reducing the cutoff preserves the pattern. For a 
ring of size 8 and a cutoff of 5 (10 for the local mode), 
the state was again dominated (Prob. 0.992) by the local 
mode, with the next contributor an excited local mode. 
Going to yet larger rings (size 10, and lower cutoff) pre- 
serves the pattern, with the first non-local mode phonon 
contributing at probability level 10 -4 . Other breather 
states yield the same result: accuracy of perturbation 
theory, insensitivity to cutoff change, and eigenfunction 
domination by local modes. 
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FIG. 2: Survival as a function of time; the initial state is the 
singly excited local mode. 

For nearest-neighbor nonlinearity the story is the same. 
Systematic study of a size-8 ring with cutoffs of 6, 8 and 
12, showed that the local mode dominated the Hamilto- 
nian's ground state. 

The message of this dominance is that the eigenstates 
of the Hamiltonian are well-localized, which is to say the 
quantized breather states are stable. 

Quantum time dependence. Although we have shown 
the breather to be dominated by the local mode, the 
point has been raised that by Fermi's Golden Rule any 
coupling to a continuum implies instability. This is not 
true. Some years ago it was found that quantum sys- 
tems with ample continuum coupling nevertheless could 
survive indefinitely |l4|. This departure from ordinary 
decay depended on smoothness and threshold properties 
of the coupling. A related matter concerns normalizablc 



bound states in the continuum |15j . 

In Fig. [3 we show \{ipo \ exp(—iHt/h)\ipo)\ 2 for tpo a lo- 
cal mode phonon and H the full Hamiltonian. As in , 
initial decay is followed by stabilization bounded away 
from zero. In [14| the amplitude was typically a few tens 
of percent, but here, consistent with Table U the ampli- 
tude stays near one. Despite the many connecting states, 
after a short time amplitude ceases to leak. The expla- 
nation is that the true eigenstate of the Hamiltonian 
has order unity overlap with the initial state. 

Note that stability may not persist in all dimensions; 
certainly threshold features of the density of states and 
spectrum are affected by dimension, and the usual in- 
tuitions regarding Fermi's Golden rule may again hold 
sway. In |8j we made the point that the symmetry break- 
ing of the Jahn- Teller effect makes this a one-dimensional 
problem, significantly enhancing the possibility of classi- 
cal breathers. The same is likely to be true quantum 
mechanically. It may even be that this plays a role in the 
temperature-dependent decay of the breather (through 
an effective increase in dimension), as evidenced by the 
high-temperature disappearance of anomalous decay in 
doped alkali halides [lj|. This issue should be address- 
able using the path integral method. 

In conclusion, we have shown that although quantum 
tunneling does in principle convert the classical breather 
into a Bloch state, when the ability to tunnel is removed, 
an initially localized pulse of energy can be trapped in- 
defnitely as a quantum excitation. 
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